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We propose a collinear factorization formula for the associated production of one 
particle and a Drell-Yan pair in hadronic collisions. It is shown that additional 
collinear singularities appearing in the next-to-leading order calculations that can 
O ■ not be factorized into parton and fragmentation functions are systematically renor- 

Q , malized by introducing fracture functions. Next-to- leading order coefficient functions 

for cross-sections double differential in the fractional energy of the identified hadron 
' and lepton pair invariant mass are presented. 
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[>. 1 I. INTRODUCTION 

o ' 
in 

The description of particle production in hadronic collisions is interesting and challenging 
in many aspects. Perturbation theory can be applied whenever a sufficiently hard scale 
characterizes the scattering process. The comparison of early LHC charged particle spectra 
with next-to-leading order perturbative QCD predictions l| shows that the theory offers 
5h 1 a rather good description of data at sufficiently high hadronic transverse momentum, 



of the order of a few GeV. For inelastic scattering processes at even lower transverse 
momentum, the theoretical description in terms of perturbative QCD breaks down since 
both the coupling and partonic matrix elements diverge as the transverse momenta of 
final state parton vanish. In this paper we will study the semi-inclusive version of the 
Drell-Yan process, Hi + H 2 — > H + 7* + X, in which one particle is tagged in the final 
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state together with the Drell-Yan pair. In such a process the high invariant mass of the 
lepton pair, Q 2 , constitutes the perturbative trigger which guarantees the applicability of 
perturbative QCD. The detected hadron H could then be used, without any phase space 
restriction, as a local probe to investigate particle production mechanisms. The evaluation 
of 0(a. s ) corrections shows that there exists a class of collinear singularities escaping the 
usual renormalization procedure which amounts to reabsorb collinear divergences into 
a redefinition of bare parton and fragmentation functions. Such singularities are likely 
to appear in every fixed order calculation in the same kinematical limits spoiling the 
convergence of the perturbative series. We therefore show how to improve the theoretical 
description providing a generalized procedure for the factorization of such additional 
collinear singularities. Most important, the latter is the same as the one proposed in 
Deep Inelastic Scattering [2]] where the same collinear singularities pattern is also found, 
confirming the universality of the collinear radiation between different hard processes. 
The latter will make use the concept of fracture functions and the renormalization 
group equations associated with them In a pure parton model approach, these 

non-perturbative distributions effectively describe the hadronization of the spectators 
system in hadron- induced reactions. We will demonstrate that the transverse-momentum 
integrated cross-section is finite and valid for all transverse momentum of the detected 
hadron, without any restriction imposed by the singular behaviour of matrix elements. 
We further note that this process is the single-particle counterpart of electroweak-boson 



plus jets associated production J4], presently calculated at nex-to- leading order accuracy 
with up to three jets in the final state One virtue of jet requirement is that it indeed 
avoids the introduction of fragmentation functions to model the final state, which are 
instead one of the basic ingredients entering our formalism. At variance with our case, 
however, jet reconstructions at very low transverse momentum starts to be challenging 6| 
and it makes difficult the study of this interesting portion of the produced particle spectrum. 



r-0 



In Ref. |7[ a first attempt was made to study the associated radiation in Drell-Yan type 



process with the aid of the concept of extended fracture functions 



The collinear factor- 



ization formula was then studied in Ref. Q where we were mainly concerned with the study 
of hard diffractive processes in hadronic collisions. The latter are a special case of associated 
production in which the tagged hadron is a proton at very low transverse momentum and 
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almost the energy of the incoming proton. In this paper we consider the production of 
a generic hadron and report all the coefficient functions for the cross-sections differen- 
tial in the invariant mass of the lepton pair and fractional momentum of the detected hadron. 

The outline of this paper is as follows. In Section [TT] we set the notation and briefly review 
the factorization of collinear singularities in the inclusive Drell-Yan process. In Section II I II 
we evaluate the 0(a s ) contributions to the cross-section in which the hadron is produced by 
the fragmentation of a final state parton. In Section IIVI we define the parton model cross- 
section for the associated production of one particle and a Drell-Yan pair in term of fracture 
functions and evaluate the corresponding 0(a s ) corrections. In Section |V] we present the 
finite cross sections in terms of renormalized fracture functions, distribution functions and 
fragmentation functions. The paper closes with a summary and some conclusions. Technical 
details and explicit formulas are collected in the appendix. 



II. INCLUSIVE DRELL-YAN PROCESS 



In this section we briefly review the factorization of collinear singularities for the inclusive 
Drell-Yan process, which will be then properly generalized when dealing with the associated 
production case. Consider therefore the collision of two hadrons H\ and H 2 of momenta P\ 
and P 2 , respectively: 

H 1 (P 1 ) + H 2 (P 2 )^ 1 *(q)+X, (1) 

where 7* stands for the virtual photon of invariant mass Q 2 ^> A.q CD and X for the un- 
observed part of the final state. The dilepton pair detected in the final state is the decay 
product of a virtual photon 7* created, to lowest order, by the annihilation of a quark and 
an antiquark with momenta p\ = X\P\ and p 2 = x 2 P 2 , respectively, and both assumed to 
be collinear to their parent hadron. The corresponding differential cross-section threfore is 
given by [n]] 



da{r) f 1 dxi f 1 dx9 <s~-^ 2 

3=1 1 



[jf l W/?W + (!««f(l-'»). (2) 



The puntiform cross-section is denoted by o"o = 4ira 2 m /9ShQ 2 ■ The hadronic centre of mass 
energy is denoted by Sh = (Pi + P2) 2 and the relevant ratio r is defined by r = Q 2 / Sh- 
To zeroth order in a s , the partonic sub-process squared energy s = (pi + p 2 ) 2 is fixed to 



i 



be s = X\X 2 Sh = Q 2 , from which the constraint w = t/xiX2 = 1 follows. The flavour 
sum runs on quarks only. The puntiform cross-sections cr is weighted by the convolution 
over parton distributions, j] 1 ' and fjj*\ evaluated at fractional momenta x\ and x% . The 
integration limits are fixed by momentum conservation. We further assign to each parton 
distribution function an index n = [1,2] denoting the hadron from which the quark (or 
antiquark) has been extracted. Calculations are performed in dimensional regularization 



with space-time dimension set to n = 4 — 2e. Following Ref. 

partonic sub-process q(pi) + q{p 2 ) - > 7*(?) is given by 

da g q(w, Q 2 ^ 



11], the cross-section for the 



5(1 



(3) 



It defines the normalization of the partonic Drell Yan cross section and corresponds to 
multiplying all matrix elements by a factor 2N C /(1 — e), iV c being the number of colours. 
The latter factor is already absorbed in the puntiform cross-section o"o in parton model 
formula, eq. (J2J). The evaluation on next-to-leading corrections is performed by computing 
the relavant real emissions and virtual diagrams. We refer to reader to Ref. ll| for the 
explicit expressions of matrix elements and further calculational details. Both real and 
virtual terms develop poles in e~ 2 which mutually cancel when these contributions are added 
so that one finally obtains 
da(r 



dQ* 



0o 



1 dxi f 1 dx 2 s—y 

X\ I -r 



x 2 



(4) 



f q 1] (x l )ff ] (x 2 ) + (q^q)\ 5(1 -w) 



+ 



2 a 
e2-K 



7T= C &( W ) 
Ztx 



In the previous equation collinear singularities appear as poles in e 1 multiplying the leading 



order 



12j splitting functions Pij(w). The adimensional factor cq appearing in eq. (j4j) reads 

« = (*&\%Z± (5) 



, Q 2 ; r(l-2e) 

where p? R indicates the renormalization scale. The subtraction of singular terms in the 
partonic cross-sections is performed by absorbing the collinear divergences into bare parton 
distributions, which in the MS scheme amounts to the following redefinition: 



" l dv r 

m = I — M(i 



u 



1 a s T(l — e) {infi 



u 



e27rr(l 



2e) V 



/4 



Pij(u) 



X 



U 



(6) 
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The renormalized distributions do depend on the scale, /if,, at which the factorization is 
performed, and their variation with respect to it gives governed by DGLAP evolution equa- 
tions 12l | . Inserting eq. ([6]) into eq. (jl]) one can explicitely check that collinear singularities, 



proportional to poles in e, do cancel. The finite result reads: 
da(r) 



0o 



1 dxi f 1 dx2 
xi Jj_ x 2 ^ 



(7) 



2vr 



C qg (w,» F /Q 2 ) 



The infrared- finite coefficient functions Cij(w) and CVj(tu), defined in the MS scheme, are 
reported in appendix [A] In phenomenological applications it is costumary to set fi F = Q 2 in 
eq. ((7]). This choice removes large logarithms of the ratio fi F /Q 2 appearing in the coefficient 
functions. At the same time, the evaluation of parton distributions at a scale fj, F = Q 2 
accounts for the resummation of such logarithms via DGLAP evolution equations. 



III. ASSOCIATED PRODUCTION 



In the present and next section we consider the next-to-simple generalization of reaction in 
eq. (jrj, namely the associated production of a Drell-Yan pair and an identified hadron H of 
momentum h: 

H^P-l) + H 2 {P 2 ) 7 *(g) + H(h) +X. (8) 

In particular we will evaluate 0(a s ) corrections to the cross-section double differential in 
the invariant mass Q 2 of the lepton pair and the fractional energy of the identified hadron 
H, defined in analogy with e + e~ annihilation process: 

_ 2h ■ (P 1 + P 2 ) _ 2Eh 

z (P + P2) 2 Vs~ h ' [) 

The last equality in the previous equations holds in the hadronic centre of mass frame. In 
this frame z is the proportional to the detected hadron energy, Eh, scaled down by the 
beam energy, y/Sh/2. In this section we consider the 0(a s ) production mechanism in which 
the observed hadron H is given by the fragmentation of a, real, final state parton and 
address such contribution as central. We label the momenta in the partonic sub-process as 
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i(pi) + j(P2) ~ > + 7*(?)j where k and g are the four-momenta of the outgoing parton 
and virtual photon, respectively. In general such a correction is expected in the form jl3| 

da c {T,z,Q 2 ) ^ f dxt f dx 2 f dp > nfl , , da ij ^' 

dQ*dz *2lj—j—j jfiMfiMD, (z/p)-^ , (10) 

ijl 

where p = 2Ei/yfSh is the fractional energy of the outgoing parton evaluated in the hadronic 
center of mass frame and it is the partonic analogue of the hadronic z. The sum runs 
over all possible partonic sub-processes and the partonic cross-sections are indicated with 
da l i~* l "t . With respect to eq. (j2J), eq. (fTUl) does contain an additional convolution on the 
fractional energy p of the final state parton weighted by the fragmentation functions Df . 
The latter gives the probability that a parton with fractional energy p fragments into the 
observed hadron H with fractional energy z. One of the diagrams contributing to eq. ffTOj) 
is depicted in Fig. (pQ). We found useful to rewrite the convolution formula as a function 
of y = (1 — cos#)/2, where 9 is its angle between the parton I and the hadron Hi in the 
hadronic center of mass frame. Within these definitions the parton-level invariants u and t 
in the matrix elements can be rewritten in this frame as 

i = {pi-kf = -Hi-w) — -p\ — - (ii) 

xiy + x 2 {l - y) 

U = [P2 - k) = - S (1 - W) r , 

x 1 y + x 2 {l -y) 

where s = (pi + p 2 ) 2 and w = r/(xix 2 ). The phase space reads: 

dp s {2) = ^{^y - T f~ 2e fay + x ^ - yrf e ~ 2 fvfr^ d y y~ e ^ - ^ e ■ ( 12 ) 
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The angular variable y and the fractional energy p of the emitted parton are not independent 
and constrained by: 

xiy + x 2 (l - y) 

The available phase space must take into account that there should be enough energy for the 
production both of the hadron H and the virtual photon 7*. This phase space constraints 
will appear in the convolution limits in eq. ffTOj) . In order to obtain them, we notice that 
the parent parton of the observed hadron is required to have a fractional energy p > z > 0. 
Applying this last constraint to eq. f TT3]) . one is able to determine the boundaries r\ and r 2 
in the X\ and x 2 convolutions integrals. They are both z and y dependent and read 

t + z(l - y) 



ri(r,z;y) 
r 2 (r,z;y) -- 



1 - zy 
r + X\zy 



xi - z(l - y) ' 

The 0(a s ) corrections in the central region therefore reads 



da c (r) 
dQ 2 dz 



, f 1 dxi f 1 dx 2 

dy — / — > 



X 2 



(14) 



(15) 



1 a s c 



fW( Xl ) f?(x 2 )D«(z/p) + (q^q) 



p2vr e" 
1 a s c 
p2n e 



P qq (w)[6(y) + 6(l-y)} 



2^ 



f?M A 1] (x 1 )Df(z/p) + (q^q) 



P q9 HS(y) + ^K qg {y) 

l« s Co- 



P ? »5(l - y) + ^-K qg {\ - y) 



The 



where p = p(y) via eq. (fl3|) and we used the shorthand Kij(y) for Kij(y,xi,x 2 ,w) 
expressions for infrared finite coefficients ifjj are collected in appendix Eq. (1151) displays 
two disjoint singular limits for y — y and ?/ — > 1. In order to expose the collinear singularites 
we have performed an e-expansion on the angular variable y: 



y) 



y 



-l-e 



-1-e 



— 5(l-y) 

€ 

--J(y) + 



1 



y 



[0,1] 

\ny 



/ ln(l-y) 
A l-y 



-[0,1] 



(9(e 2 



2/ '+[0,1] 



G(e 5 



+[s,i] 

The unregularized splitting functions appearing in eq. (TTBl) are given by 

1 + w 2 



T R [w< 



w 



(16) 
(17) 

(18) 



with Cp = 4/3 and Tr = 1/2. We wish to conclude this section by noting that the collinear 
divergences appearing in eq. f lT5"j) do correspond to configurations in which the parent parton 



s 



of the observed hadron is collinear to the incoming parton. Such divergences at vaninish- 
ing transverse momentum escape, as shown already in the context of Semi-Inclusive Deep 



Inelastic Scattering {2], any factorization in terms of renormalized parton distributions and 
fragmentation functions. While for many practical applications they are regularized intro- 
ducing an arbitrary cut-off on the produced hadron transverse momentum, configurations 
which give rise to these divergences will be present at every order in perturbative calcu- 
lations. Fracture functions together with their own renormalization group equations can 
be shown to provide the correct tool to perform the resummation to all orders of large 
logarithmic contributions coming from the factorization of such collinear singularities. 

IV. CORRECTIONS IN THE TARGET REGION 

To lowest order in the QCD coupling no hadron can be produced in the final state since 
QCD radiation is absent. In this case we assume that hadron production is described by 
fracture functions (xi,z) and M^ H2 (x2, z). These non-perturbative distributions 

give the conditional probability of finding a parton of fractional momemntum X\(x2) in 
the incoming hadron Hi(H 2 ) while an hadron H, with fractional momentum z, is detected 
in the final state jjj. In a pure parton model approach they describe hadron production 
in the target fragmentation region of H\ or H2. The latter regions, denoted by TZt ± and 
1Zt 2 i respectively, can be defined as 9 = and 6 = tt, where 9 is the angle between H and 
Hi defined in the centre of mass frame. Fracture function were originally introduced to 
describe hadron production in the Deep Inelastic Scattering target fragmentation region. A 
renormalization group evolution equations were derivered [31 with the aid of Jet Calculus 

n 

technique 114[. Subsequently the soft and collinear factorization of these distributions in 
Semi- Inclusive DIS wasproven respectively in Refs. j3> Q] • A complete one loop calculation 
was presented in Ref. {2], confirming the factorization conjecture first formulated in Ref. {jj. 
We emphasize that by virtue of the factorization theorem, fracture functions are univeral 
distributions, at least in the context of Semi-Inclusive DIS. In hadronic collisions however 
such a proof does not exist and_counter example to it have been given in the context of 
diffractive production in Ref. 



15| . The main motivation for using fracture functions in 
the present context is that, as we shall prove in the following sections, they allow us to 
systematically factorize collinear divergences occuring in the evaluation of the partonic cross- 
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FIG. 2: Pictorial representation of the parton model formula in eq. (|19p . The hadron H is non- 
perturbatively produced by a fracture functions. 



sections. Given these assumptions, we now present the parton-model formula, depicted in 
Fig. ([2]), for the associated production case: 

da H (r) 



dQ 2 dz 

z dx\ f 1 dx2 



(19) 



f~ Z — [ ^^ej[MW(x 1 ,z)/f ] (^) + M[ 11 (x 1 ,z)4 2 ](x 2 )l5(l-t i ;) + 

J T Xi J^r_ X2 ~~ L J 



J —L — 1 J -2L- $2, 

i—2 xi q 



The integration limits of convolution integrals in both lines of eq. f TT9]) are given by mo- 
mentum conservation: 

1 — z > X\X 2 > t and X\ + z < 1 in IZ^ , (20) 

1 — z > X\X2 > t and X2 + z < 1 in 1Zt 2 ■ (21) 

Phase space integrations are asymmetric since each fracture function selects its own frag- 
mentation region. To uniform the notation, we exchange the superscript H/H\ (H/H2) for 
[1] ([2]) which proves to be useful for the bookkeeping of the various contributions. Fracture 
functions in eq. (fT9~l) are normalized according to the constraint 

V / dz daH{T) - ^ (22) 
^ I dZ dQ 2 dz ~ dQ* ■ 



H 



The above constraint must me fullfilled irrespective to the order of the perturbative calcu- 
lations. It is interesting to note that for the inclusive Drell-Yan cross-sections appearing on 
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the right hand side of eq. (122 p . the factorization theorem guarantees that the corresponding 
cross-sections can be described by universal parton distributions functions. On the left hand 
side of eq. (1221) . we have instead no guarantee that fracture functions eventually extracted 
from Deep Inelastic Scattering can be succesfully used in hadronic collisions. Eq. ( TT9|) is 
therefore both a factorization and a cross-section conjecture for the process under exami- 
nation. In the remainder of this section we consider the evaluation of 0(a s ) corrections to 
the parton model formula, eq. fTT^|) . We address it as target contributions to distinguish 
them from the one evaluated in Sec. IIHI As already stated, when the final state hadron is 
observed in IZ^ or TZt 2 , we assume that it has been non-perturbatively produced from a 
fracture functions. Final state partons occurring in 0(a s ) corrections to eq. ffT9l) must be 
therefore integrated over phase space and virtual corrections added. A diagram contributing 
to eq. f[2~3"l) is depicted in Fig. ([3]). The calculation closely follows the one already presented 
in Sec. [TT] for the inclusive Drell-Yan process upon the exchange of a parton distributions 
with a fracture functions and taking into account momentum conservation in integrations 
limits. We therefore just quote the final result, valid up to 0(a s ): 

da T (r) f 1 ^ z dxi f 1 dx 2 



dQ 2 dz 



0o 




M q l \x l)Z )ff\x 2 ) + (q O q)\ [5(1 - w) 
(Mjpfa, z)ff\x 2 ) + M q l \x 1 , z)ff\x 2 )) + (q^q) 

1 a 



(23) 



" S P U W ) c + 7T-C qq -{w) 



€ 2tt 



e 2tt 



Pqgiw) C 



2tt 



2^ 



C qg{ w ) 



+cr 



dx\ 



T X\ J J_ 



1-2 



dx-2 
x 2 



Mf\x 2 , J)/! 1 ^,) + (q « g)] \S(1 - w) - ~P„(w) c„ 

H J L € ZTT 

{Mf\x 2 , z)f q 1] (x 1 ) + Mf\x 2 , + (q qj 

1 a 



Z7T 



+ 



a. 



-^-Pqq{ W ) C + 7^-C< 



6 27T 



2tt 



19 



W 



Comparing eq. (jlj) with eq. (123]) reveals that in target fragmentation region the structure of 
collinear singularities is the same as in the inclusive Drell-Yan expected. The main 

change is just a restriction on phase space integrals since the production of the Drell-Yan 
pair of a given invariant mass Q 2 must be, by energy-momentum conservation, compatible 
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FIG. 3: Example of 0(a s ) real corrections contributing to eq. (|23p . The hadron H is non pertur- 
batively produced by a fracture function. 

with the observation of a hadron in the final state with fractional momentum z. 

V. FINITE CROSS-SECTIONS AT NLO 



In this section we will describe the collinear factorization procedure which must be applied in 
order to get infrared finite results for the cross-section under examination. We have already 
noted, by comparing eq. (J4J) and eq. (|2"51) . that the structure of collinear singularities in the 
target fragmentation region is identical to the one found in the inclusive Drell-Yan case. 
We may expect that renormalized fracture functions are defined in a way similar to that of 
renormalized parton densities in eq. ([6]). As it was firstly obtained in the original analysis 
of Ref. {2! and confirmed in the one loop calculation of Ref. j^j], the renormalized fracture 
functions obey a somewhat more involved subtraction with respect to parton distributions. 
In the MS scheme the redefinition of bare fracture functions reads: 



Mf IH \x,z) 



du 



x/(l-z) 



u 



8ij8(l — u) + - 



la, rfi 



e27rr(l - 2e] 



/4 



>«(«) 



M, 



H/H 



+ 



x/{x+z) du 1 ula s T(l - e) / 47174 



u 1 — u x e 2n T(l — 



2e) V 



p ij( u ) fj/HA-) D 



X 



u 



zu 



X 1 



u 



(24) 



where in our notation Pij{u) = P(i)i^.j{u). The first term on r.h.s of eq. f[24"|) has the same 
subtraction structure as for parton distribution, eq. ()6]). The singularity is due to collinear 
radiation accompaining the active parton, while the hadron in the final state is non perturba- 
tively produced by the fracture functions itself. In the second term of eq. (124]) the singularity 
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is due configurations in which the parent parton of the observed hadron is collinear to the 
incoming parton. The factorization procedure is accounted for by substituting in eq. ( j2SJ) the 
bare fracture and distributions functions by their renormalized version in eq. ([6]) and eq. (J24]) . 
Renormalized parton distributions and fracture functions homogeneous terms do cancel all 
singularities present in eq. fl23|) . The additional singularities in eq. (1151) are cancelled by the 
combination of parton distributions and fracture functions inhomogeneous renormalization 
terms. The final result, up to order 0(a s ), is obtained adding the the various contributions: 



da M r 



da 



T 



da H > c \ 



T 



dQ 2 dz dQ 2 dz 
l ~ z dxi f 1 dx 2 



dQ 2 dz 



00 



E 



(25) 



MW(x l5 z^l)ff\x 2 ^ 2 F ) + (q^ q)] [s(l - w) + ^C«(w, ^) 



+ 



1 dxi f x ~ z dx 



O.s 



/4 



T X\ J J_ 

1 — 2 



X-2 



Mf\x 2 ,Z,li 2 F )f^{x u iil) + {q^q) 5(1 -w) + ^C q ,(w 



'' Q 2 



+ 



(Mf\x 2) z^l)ff(x l ^l) + Mf\x 2 ,z^l)ff\x 1 ^ 2 F )) + (q O q) 



J 2tt 



/4 



— C qg [W, — 



1 d y ^^1 
Jri 



'■_> 



x 2 



4 1] (*i,/4) /fWl) + (q » q)] Df(z/p,^K gq {y, ^) + 



fS- ] {*u&) [fl%2,» F M(z/p,n F ) q)] ^K qg (y, ^) 



+ 



4 2l (*2,/4) fP(x 1 ,p 2 F )Df(z/p,p 2 F ) + (q^q) 



/4 
Q 2 



We have used the shorthand Kij(y, fi 2 F /Q 2 ) for Kij(y, xi,x 2 ,w, [i 2 F /Q 2 ). The explicit form of 
the finite coefficient functions K qq and K qg is reported in appendix^] In the last three lines of 
eq. ( 125]) we let depend parton distributions and fragmentation functions on the factorization 
scale since this replacement induces subleading corrections to the current accuracy. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have calculated the 0(a s ) corrections to the associated production of one 
particle and a Drell-Yan pair. Additional collinear singularities found in the perturbative 
calculations do correspond to configurations in which the parent parton of the observed 
hadron is collinear to the incoming parton. These singularities can, to 0(a s ), be consistently 
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absorbed into renormalized fracture functions and resummed to all orders by using the 
evolution equation given in Refs. as. With this technique, the presented cross-sections 
does not require any cut in the transverse momentum of the observed particle, while a 
perturbative treatment is guarantee by the presence of high invariant mass dilepton pair. 
Quite imprtantly, the factorization of collinear singularities in the present context makes 
use of the collinear subtraction structure already defined in the context of Deep Inelastic 
Scattering. Despite the fact the the full results make use of fracture functions and threfore a 
phenomenological modelling of the latter would be eventually required, the advatages reside 
in that fracture functions embodies the correct scale dependence through their own evolution 
equations. 



Appendix A: Finite coefficients 



In this Appendix we present the results for the finite coefficient which appear in the 
previous sections. The plus distribution are defined in the usual way: 



1 — w 



(Al) 



'o (1 - W)+[ ,1] " JO 

The subtraction point is underlined. The coefficient functions in the target fragmentation 
region do coincide with the one found in the inclusive case and read: 



C qg [w, 



Q 2 J 



-2P qq (w)ln^ + C qq (w) 
-P qg ( w )ln^ + C qg (w), 



(A2) 



where the scale independent coefficients Cy are given by: 

— / Inf1 — ml \ 

C qq (w) = C F 



1 — w 



-[o,i] 



2 1±^! In w + (-7T 2 - 8)5(1 -w) 
1 — w V3 



C qg {w) = T R 



(w 2 + (1 - wf) In 



:i-w) 



w 



1 



7 



- + 3w — —w' 



(A3) 



The polynomial term in C qg is slightly different from the one reported in Ref. 11] because 
an additional term (1 — e) _1 is provided for matrix elements with a gluon in the initial 
state. This accounts for the correct gluon polarization in n-dimensions and it affects only 
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the polynomial terms in the coefficient function. For the central term the subtraction is 
performed on the angular variable y. The plus distributions in eq. (fTtjj) are defined by: 



V h(v) 



o (1 - 2/)+[o,i] Jo 1 - y 



1 . My) 

"2/7! = 

lyj+[o,i] Jo 



(A4) 
(A5) 



where is a smooth test function. The coefficients Ky read: 
Q 2 



K q q[y,x 1 ,x 2 ,T,^)=K qq -(y,x 1 ,x 2 ,T) + 



-<%)Pg 9 0) — In & - 5(1 - y)P qq (w)—^\n * 



XiX 2 — r Q 1 



x x x 2 - T Q 2 ' 



/4\ 



2 

( y, an, ^2, r, ^ J =K qg (y, x lt x 2 , r) - <%) ^— P g9 (w) In . (A6) 



Q 2 



By defining b(xx,x 2 ,y) = x±y + x 2 (l — y) and a = x\x 2 — r, the scale-independent 
coefficients read: 

t> , \ „ \ l ( X2 r, \ 2rb2 \ 1 1 f x i 2rb2 
K qq {y, x u x 2 , r )=C F < - — (1 - y) + — + t V + -77 , /, x 

\b\Xi a 2 / 2/+[o,i] b \ x 2 a J (1 - 2/)+[o,i] 



+ S(y) 
K qg (y,x 1 ,x 2 ,r)=T R 



1 1 + w 2 x 2 . x\t 

In — 

x\ 1 — w a er 

x 2 2(1 -y)A 1 



+ S(l-y) 



ay 



1 1 + W 2 X\ . x\t 

— 1 ln — 

x 2 1 — w a a z 



bx\ J y+\o,i] x 2 b 2 



x 2 



w 2 + (1 - wf 



111 



TXn 



1 



(A7) 



The coefficient function K gq can be obtained by exchanging x\ -B- x 2 and y -H- 1 — y in 
the expressions for In one can note the appearance of poles in (1 — w)' 1 , where 
w = t/(x\x 2 ). In this limit only soft parton emissions are allowed. This limit however is 
outside the integration region specified in the central contribution to the final result, eq. (|25|) . 
due to the requirement z > embodied in the specific form of the integration boundaries r\ 
and r 2 in eqs. (TT4T) . 
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